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I. INTRODUCTION 


The optical theorem is a basic result in scattering theory that is of both fundamental 
interest and considerable applied importance [T]. It can be formulated in a variety of settings, 
including quantum mechanics, acoustics and electromagnetic theory. In its simplest form, 
the optical theorem relates the power extinguished from a plane wave incident on a scattering 
medium to the scattering amplitude in the forward direction of scattering. For the case of 
electromagnetic scattering, the extinguished power Pe is given by 

Pe = -^ImA ■ E* , (1) 

2fco ° ^ ^ 

where ko is the wavenumber, Eq is the incident field and A is the scattering amplitude in 
the forward direction M. In physical terms, the loss of power from the incident field is 
due to interference between the incident held and the scattered held within the volume of 
the scatterer. 

The standard derivation of the optical theorem makes use of the ansatz that the scattered 
held Es behaves asymptotically as an outgoing spherical wave of the form 

^ikor 

Es ~ A- , kor —>-00. (2) 

r 

The above ansatz may be justihed for the case of material media in which the polarization is 
related to the electric held by a linear constitutive relation. Moreover, in this situation, the 
optical theorem may be derived without invoking the asymptotic behavior of the scattered 
held [5l E] . 

The optical theorem is normally considered within the framework of linear optics |1]. 
However, the ansatz (|^ is very general. That is, all properties of the scatterer are encoded 
in its scattering amplitude, which can, in principle, be arbitrarily prescribed. Thus, there 
is no reason to restrict the optical theorem to the linear response regime. Making use of 
this observation, in this paper we consider the optical theorem in the context of nonlinear 
media. We derive an expression for the extinguished power that holds when the polarization 
is an arbitrary function of the electric held. To some extent, this result may be expected 
on physical grounds. However, we provide a proper mathematical justihcation following the 
approach of nu- We specialize our result to the cases of quadratic and cubic nonlinear- 
ities. We also study in detail the processes of second harmonic generation and the Kerr 
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effect for small scatterers. Our results on scattering from small nonlinear particles are of 
independent interest, since exact solutions to nonlinear scattering problems are, to the best 
of our knowledge, known only in one dimension 13 E]. 

The remainder of this paper is organized as follows. In section II, we present some basic 
results in nonlinear optics and scattering theory that will be used later in the paper. Section 
III presents the derivation of the optical theorem in the form we require, without the use 
of asymptotics. In section IV, we consider separately the cases of second- and third-order 
nonlinearities. Numerical results for small scatterers are presented in section V. Finally, in 
section VI we consider applications to apertureless scanning near-held optical microscopy 
in which the tip exhibits a nonlinear optical response. We investigate the cases of second- 
and third-harmonic generation and characterize the achievable resolution for a model system 
consisting of two scatterers. A discussion of our results is presented in Section VII. 

II. PRELIMINARIES 

In this section we collect several results in nonlinear optics and scattering theory that 
will be useful in the derivation of the generalized optical theorem. We begin by recalling 
that the Maxwell equations in a source-free nonmagnetic medium are of the form 


(3) 

(4) 

(5) 

( 6 ) 


V-D = 0 , 



V-B = 0 , 



Here E is the electric held, B is the magnetic held, D is the electric displacement held and 
P is the polarization. In addition, D and P satisfy the relation 


(7) 


D = E -|- dvrP . 


Throughout this paper, we will use the following Fourier transformation convention: 



( 8 ) 


( 9 ) 
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where the time and frequency dependence of all quantities are displayed explicitly. We note 
that if /(r,f) is real-valued, then /(r, —cj) = f*{ir,u). Upon Fourier transforming (|^-([^ 
we obtain 



<1 

U 

0 

(10) 

V X E(r,a;) 

— ik{u)'B{r, u) = 0 , 

(11) 


V • B(r,ci;) = 0 , 

(12) 

V X B(r,C(;) 

+ ik{uj)Y){r, cj) = 0 , 

(13) 


where k{co) = oj/c. 

If the medium is linear, the polarization is given by 

Pi{r,u}) = Xif{r-,u)Ej{r,u) , (14) 

where Xi]'^ is the hrst-order susceptibility. Here we have adopted the summation convention, 
whereby repeated indices are summed. For the case of quadratic nonlinear media. 


= x[]\r;uj)Ej{uj) + ^ xJ2(r;Wi,a;2)Uj(r,n;i)Ufc(r,^ 2 ) , (15) 

( 2 ) 

where Xi/k i® second-order susceptibility. The sum implies that the electric held at the 
frequencies uji and UJ2 contributes to the polarization at the frequency oj ii uji+ UJ2 = oj- For 
cubic nonlinear media. 


Pi{r,uj) = xlf{r]u)Ej{u) + ^ xfjliir]UJi,u}2,u}3)Ej{r,Ui)Ek{r,U2)Ei{r,U3) (16) 

UJ-i+UJ2+l^3=‘^ 

where xifkii'^'i ^ 2 , ^ 3 ) are the third order susceptibilities. 

We will assume that the susceptibilities have over all permutation symmetry. This as¬ 
sumption is quite standard and holds for nonresonance frequencies in the classical anhar- 
monic oscillator model and in quantum optics Ei. However, we note that Kleinman 
symmetry can be broken in a variety of systems BM- 

The wave equation for the electric held E(r, u) is obtained by taking the curl of (|^ and 
eliminating the magnetic held B(r,a;) using (§. We then have 

V X V X E(r,n;) - k‘^{uj)E{r,uj) = 47rk‘^{uj)P{r,u:) . (17) 


The solution of the wave equation 0 is given by 

Ei{r,uj) = Einc,iir,u) + k‘^{u) / £rGij{r,r'-,u)Pj{r',u), 


(18) 
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where i^) obeys (17) with P = 0. The Green’s function Gij is of the form 

Gij{r,r';uj) = G{T,r';u) , (19) 


where 


G(r,r';a;) = 


gifc(a;)|r-r'| 


r — r 


( 20 ) 


Straightforward calculation shows that the Green’s function has the following asymptotic 


form 


Gij{r,r';u) 


^ik{Lj)r 


A 


t] Tirjje 


—ik{Lj)r-r' 


( 21 ) 


when r ^ r' . Using this result, we hnd that the scattered held behaves as an outgoing 
spherical wave of the form 


^ik{u])r 

P^,i(r,a;)- Aiir.u) , 

r 

where the scattering amplitude is dehned by 

Ai{r,bj) ^ - TiTj) [ 


( 22 ) 


(23) 


Following standard procedures, the conservation of energy follows immediately from (17) 
and takes the form 


V • S(r,a;) = ^^^^Im(E*(r,w) • P(r,a;)) , 
where the Poynting vector S is dehned by 

= ^Re(E(r,a;) x B*(r,a;)) . 

We recall that the time-dependent Poynting vector is dehned as 

S(r,t) = ^E(r,f) X B(r,f) . 
dvr 

For time-harmonic helds, the time average of the Poynting vector 


is well dehned. It follows that 


S= hm;^/' S{r,t)dt 

T^oo 1 Jq 


S = / S{r, u)du . 
Jo 


(24) 


(25) 


(26) 


(27) 


(28) 


We note that for nondispersive media, the analog of the Manley-Rowe relations can be 
own to hold. That is, if the susceptibiliti 
V • S = 0. The proof is given in Appendix A. 


shown to hold. That is, if the susceptibilities Xi]\ xf^k xfjli purely real, then 
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III. OPTICAL THEOREM 


In this section we derive the optical theorem for nonlinear media following the approach 
of in E]- We begin by recalling some basic facts from scattering theory |3]. We consider 
a general nonlinear medium, whose polarization is defined by either a quadratic or cubic 
nonlinearity. We suppose that a held Ej„c is incident upon the medium and write the total 
electric held as the sum 

E = Einc + E^ , (29) 

where E^ is the scattered held. It follows from ( [ItI ) that the scattered held obeys 

V X V X Es(r,ci;) — /c^(a;)Es(r,a;) = 47rk‘^{uj)P{r,u) . (30) 


The energy carried by the scattered held is associated with the Poynting vector S<j, which 
is dehned by 


S. = ^Re(E.(i..) X B:(ix)) . 


(31) 


Evidently, obeys the conservation law 


V ■ S,(r,ix) = d^Im(E;(r,ix) • P(r,x))) , 


(32) 


which is a consequence of (30). 


Suppose that the scattering medium is contained in a volume V. Then the power absorbed 
by the medium is given by 


Pa{uj) = - S(r,a;) ■ hcfr , 

JdV 


(33) 


where h is the outward unit normal to dV and the presence of the overall minus sign signihes 
that this is the hux of the Poynting vector of the wave traveling into the medium. Converting 
the above surface integral to a volume integral by means of the divergence theorem and 


making use of (24), we have 


ckioj) 


Pi.(x)) ^ J Im(E"(r',x)) ■ P(r',ix))<iV . 


(34) 


In a strictly similar manner, we dehne the scattered power as 


Ps{uj) = I Ss(r',a;) ■ h(fr' 
Jav 


(35) 
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We then obtain from (32) that 


PJoo) = 


ck{oj) 

27r2 


Im(E*(r',a;) ■ P(r',. 


(36) 


'V 


We dehne the extinguished power to be the total power lost from the incident held due 
to absorption and scattering: 

Pei^) = Pa{^) + Psi^) ■ (37) 


It follows from (34) and (36) that the extinguished power is given by 

ck(u>) 


PM = 


27r2 


'V 


■P(r',u;))dV . 


(38) 


We note that if either Eqs. (15) or (16) for the polarization is inserted into the above 


expression, we see that the power extinguished from the incident held is due to interference 
between the incident held and the total held within the volume of the scatterer. 


We can now rewrite (38) in terms of the scattering amplitude, provided that the incident 


held is a plane wave of the form 


E.„,(r,i*)) = E„(i^)e“l">‘-' 


(39) 


where s is the direction of propagation. Upon comparing (38) and (23), we obtain the 
optical theorem 


PM = 


-ImA(s, (j) ■ EQ(a;) . 


27r^k(cj) 

Using this result, we see that the time-averaged extinguished power is given by 

bOO 


/•OO ^ /•c 


k(uj\ 


-ImA(s,a;) • EQ(a;)(icu . 


(40) 


(41) 


We note that the optical theorem (40) applies to both linear and nonlinear media. In 


Section IV, we specialize this result to the cases of quadratic and cubic nonlinearities. Here 
we remark that in the case of a linear medium with an incident monochromatic held of the 
form 

E,„,(r, (u) = - H) + E*<5(u; H)) , (42) 


(41) becomes 


2M(kl) 


ImA(s,U) -E* 


(43) 


We have thus recovered the familiar form of the optical theorem ([^. Finally, we note that 


(40) is an exact result. That is, it has not been derived by making use of the asymptotic 


behavior of the electric held. 
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IV. SECOND- AND THIRD-ORDER NONLINEARITIES 


Evidently, in order to apply the optical theorem (41) it is necessary to hrst obtain the 


scattering amplitude. In this section we calculate the scattering amplitude for second- and 
third-order nonlinearities. We begin with the case of second-order nonlinearity and, for 
simplicity, discuss only the problem of second-harmonic generation (SHG). In this setting, 
we analyze the scattering of an incident monochromatic wave from a spherical particle 
whose size is small compared to the wavelength. Next, we turn our attention to the case 
of third-order nonlinearity, where we restrict our attention to the Kerr effect. Once again, 
we calculate the extinguished power for a small particle and study the associated resonant 
scattering. 

We note that the method we develop for calculating the scattering of light from a small 
nonlinear inhomogeneity may be of independent interest. In particular, it is readily extended 
to collections of small inhomogeneities, which is a physical setting that arises in applications 
to biomedical imaging and nonlinear microscopy [THHIS]- We plan to report the results of 
such calculations elsewhere. 


A. Second Order Nonlinearity 

We consider SHG excited by a monochromatic incident held of frequency G. We assume 
that the second-order susceptibility is sufficiently weak so that the condition 


X^^l{r;uJi,uJ2)Ej{T,uJi)Ek{r,uj2) < i^;^)Ej{uj) 


( 1 )/ 


(44) 


U)\^UJ2=^ 


is obeyed. We then hnd that the wave equation 0 together with (|l^ and the permutation 


( 2 ) 

symmetry of Xi/k gives rise to the pair of coupled wave equations 

V X V X E(r,G) - k\n)E{r,n) = 47rA;2(G)(yW(r, G)E,(r, G) + 2xg(r, 2G,-G)E,(r, 2G)E*(r, G)) , 


(45) 


V X V X E(r,2G) - k'\2Q)E{r,2Q) = 47rF(2G)(xy^(r, 2G)E,(r, 20) + xg^(r, G, G)E,(r, G)Efc(r, G)) , 


42), 


(46) 


for the electric helds at the frequencies O and 20. Note that we have not accounted for the 


formation of higher harmonics, consistent with the condition (44) 





It follows immediately from (18) that the solutions to (45) and (46) are given by 

Ei{r, fi) = Einc,i{T^, fl) + k^{Q) j dVxJJ.^(r', r'; fl)Efc(r', fl) 

+ 2k\n) j dV'xgi(r', 21^, -f2)Gi,(r, r'; n)E,{r', 2n)E*{v\ fl) , (47) 

E,{r,2Q) = k\2Q) j d^r'xf,!{r',2Q)G,,{r,r'-2Q)Ek{r',2Q) 

+ k\2Q) j dV'xgJ(r', ^^)G'^,(r, r ; 2Q)Ek{r', Q)Ei{r', Q) . (48) 

Suppose that the scattering medium is a small ball of radius a with k{Q)a <C 1 and 
/c(2fl)a ^ 1. The susceptibilities are taken to be = dif x^^^(r;a;) = ri^jl for 

|r| < a and to vanish for |r| > a. Eqs. (47) and (48) thus become 


(49) 



+ k\^)Vjk J 

f d^r'G,^ 

(r, r';G)T;fc(r',l 


|r^|<a 


+ 2k\Q)gf^i 


(r,r';G)Efc(i 

■',2G)E;(r',G) 

' |r^|<a 



= k‘^{2Q)gf^ 

[ dh'G.ji 

(r,r';2G)Efc( 

y,2G) 

' |r^|<a 



+ k‘^{2n)gll\ 

[ dh'G.ji 

' |r^|<a 

(r,r';2G)Efc( 

y,G)y(r',G) 


(50) 


Using the asymptotic form of the Green’s function given in (21), we hnd that the scattered 
helds are of the form 

^ik(Q)r 


Et{r,n) = A{r,n)- 

Ei{T,2^) = Ai(r,2G) 


where the scattering amplitudes are dehned by 


r 

^ik{2Q.)r 

r 


(51) 

(52) 


dvr 


A(r,G) = ^AAj-fA)k'\n)^l^Ek{0,n) + A^ZEk{0,2n)E:{0,n)), (53) 


( 2 ) 


dvr 


7l,(r, 2G) = ya^(5.,- - fA)e{2n){T]fjEjiO, 2G) + v-tiMO, ^)EiiO, G)) . (54) 

Here we have used the fact that the radius of the scatterer is small, which leads to the 
identity 


( 2 ) 


dvr 


' |r^|<a 


4V'e-“<")'-''g(l;(w)r') = —o=‘g(0)(l + 0{k{L])a)) , 


(55) 


for some function g. The local helds i7j(0, G) and i7j(0, 2f2) can be calculated perturbatively 
as shown in Appendix |B 1[ 
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B. Third Order Nonlinearity 


The treatment of the Kerr effect parallels that of SHG. We consider the Kerr effect 
excited by a monochromatic incident held of frequency O. We assume that the third-order 
susceptibility is sufficiently weak so that the condition 


X%{r;ui,U2)Ej{r,Ui)Ek{r,U2)Ei{r,U3) ^ xif {r;u)Ej{u) 


(1)/ 


(56) 


UJi+LJ2+‘^3=<^ 


( 3 ) 

is obeyed. We then hnd that (17) together with (16) and the permutation symmetry of Xijki 
gives rise to the wave equation 

V X V X E(r,0) - fc2(0)E(r,0) =ATTk‘^{n){x[f {r,n)Ej{r,n) 

+3xUi{r,n,n,-n)Ejir,n)Ek{r,n)E:{r,n)) . (57) 

Note that we have not accounted for the formation of higher harmonics, consistent with the 


condition (56). It follows immediately from (18) that the solution to (57) is given by 

Ei{r, G) = Einc,i{r, G) k^{kl) J d^r'xfki^'^ kl)Gij{r, r; kl)Ek{r', G) 

+ 3k\n) J dVxf2™(r',G,G,-G)Gi,(r,r';G)Efc(r',G)EKr',G)E;;(r',G)) ,(58) 

Suppose that the scattering medium is a small ball of radius a with k{fl)a -C 1. The 
susceptibilities are taken to be Xj'j^(r;a;) = and xifkii^'i^) = vlfki kl — ® 
vanish for |r| > a. Eq. (58) thus become 

Ei{r, G) = Ei„c,i(r, G) k‘^{Q)r]ff} 


|r'|<a 


dVGij(r, r'; G)Efc(r', G) 




|r'|<a 


dVG,,(r,r';G)Efc(r',G)Ez(r',G)E;;(r',G)) . (59) 


Using the asymptotic form of the Green’s function given in (21), we find that the scattered 
held is of the form 


Ef(r,G)=A(r,G)- 


where the scattering amplitude is dehned by 


(60) 


/l.(r,SJ) = + 3),«„Et(0,n)E,(0,n))B;(0,n)) ,(61) 


.(3) 


Once again, we calculate the local helds Ej(0, G) perturbatively, as shown in Appendix |B 2[ 
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V. NUMERICAL RESULTS 


In this section we apply the optical theorem (41) to linear, second- and third-order non¬ 


linear media. We present numerical results for several cases of interest, including second 
harmonic generation and the Kerr effect. We will see that the effect of the nonlinearities is 
to modify the linear scattering resonance of small scatterers. 


A. Linear Response 

We consider an isotropic medium with = T]^^^6ij and assume that all the higher-order 
susceptibilities vanish. The incident held is taken to be a unit-amplitude plane wave of the 
form = Eoexp(i/coS-r) with Eq = Eqx and s = z. Using (98) and (91), the extinguished 
power becomes 

SO _ „ / nW \ 

' (62) 


8U o o 


We can write the above formula in a more familiar form in terms of the renormalized polar¬ 
izability a, which is dehned as 


a = 


ao 


(63) 


1 — [l/a -|- i^k] ’ 

where ao is the zero-frequency polarizability, which is dehned in terms of the linear dielectric 

6 ( 1 ) 


permittivity e(i). Here 


1 


ao — 


41 ) + 2 


where e(i) = 1 -|- We hnd that (62) becomes 

2U 


Pf, = —Eglma . 
TT 


(64) 


(65) 


In Figure 1 we illustrate the frequency dependence of the extinguished power for a dielectric 
scatterer of size a = lOOnm with e(i) = —5.28. We see that there is a scattering resonance 
at ka ^ 0.7. 


B. Second Harmonic Generation 

We consider the case of a medium with isotropic 77 ( 1 ) and obeying permutation sym¬ 
metry. That is, and with the other rj^jl vanishing. We also assume 
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PJPq 



FIG. 1. Frequency dependence of extinguished power for a single linear scatterer. Here Pq = cl^cEq. 


that the incident held Ej„c points in the x-direction and the direction of observation s is 


taken to be in the ^-direction. It follows from (B20) that the extingnished power Pe is given 
by 


D 8^ 3t 
Pe = —a Im 
3 


{r,mEf\o,a) + ( 0 , n) + SJ))-£)P( 0 , a) 


( 66 ) 


In Fignre 2 we illnstrate the freqnency dependence of the extingnished power for SHG. 
The scatterer size is a = lOOnm and Eq is taken to have nnit amplitnde. Plots are shown 
for e = = 0, 0.1, 0.2, 0.3. We see that the resonance shifts to lower freqnencies 

and its amplitude increases by more than a factor of two relative to the linear case. The 
physical situation considered corresponds to the material / 9 -BaB 204 [8]. We note that there 
does not appear to be a simple physical argument to predict the extent or direction of the 
reported frequency shifts. 
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FIG. 2. (Color online) Frequency dependence of extinguished power for a single nonlinear scatterer 
with SHG. Here Pq = o^cEq. 

C. Kerr Effect 


We consider the case of a medium with isotropic and obeying the permutation 
symmetry. That is, and with all other riljh vanishing. We also 

assume that the incident field Ej„c points in the x-direction and the direction of observation 


s is taken to be in the z-direction. It follows from (B32) that the extinguished power Pe is 
given by 


i(0,SJ 


(67) 


In Figure 3 we illustrate the frequency dependence of the extinguished power for the 
Kerr effect. The scatterer size is a = lOOnm and Eq is taken to have unit amplitude. Plots 
are shown for e = = 0,0.01,0.02. We see that the resonance shifts to higher 

frequencies relative to the linear case. As may be expected, the effect is less pronounced 
than in the case of SHG. As above, we do not know of a simple physical argument to predict 
the extent or direction of the reported frequency shifts. 
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FIG. 3. (Color online) Frequency dependence of extinguished power for a single scatterer with 
Kerr nonlinearity. Here Pq = u'^cEq. 

VI. APPLICATION TO NEAR-FIELD MICROSCOPY 

Near-field scanning optical microscopy (NSOM) is a widely nsed experimental tool to 
overcome the diffraction limit of optical microscopy [I6]. In a typical experiment, an aper- 
tnred probe (often a metallic coated optical fiber) is bronght into the near-held of a sample 
and employed as an optical sonrce. The image is formed by scanning the probe and record¬ 
ing the intensity of light scattered into the far-held. In a reciprocal arrangement, the probe 
may be nsed as a detector with the illnmination incident from the far-held. In either case, 
the resolntion of the resnlting image is controlled by the size of the probe rather than the 
wavelength of light. 

Apertnreless NSOM is an alternative to the above approach in which the illnmination and 
detection both take place in the far-held uniE]. The experimental setup is illustrated in 
Figure]^ in which an incident held illuminates a metallic tip that is placed in the near-held 
of the sample. The image is obtained by scanning the tip and measuring the scattered held 
with a detector that is placed in the far-held of the sample and the tip. 

A rehnement of apertureless NSOM is to introduce a huorescent tip, which allows for 
the spectral isolation of the detected light and improvement in SNR by background sup- 
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FIG. 4. Illustrating the apertureless NSOM experiment. 


pression [18] . Spectral isolation may also be achieved by utilizing a tip that has a nonlinear 
optical response [I9H25|. This approach has the advantage that it is not affected by fluo¬ 
rescent photobleaching. Experiments in which SHG, THG and four-wave mixing have been 
utilized for aptertureless NSOM in a dark-held conhguration have recently been reported 


In this section we develop a model for SHG and THG apertureless NSOM. We consider 
a system in which a nonlinear metallic tip is placed in the near-held of a pair of small 
dielectric particles, which are taken to have only a linear optical response. The setup is 
shown in Figure and the mathematical details are presented in Appendix G. We begin 
with the case of SHG. The sample consists of a pair of dielectric spheres of radius A/(107r) 
and susceptibility fj^j^ = OASij, located at the positions ri = (0, 0) and r 2 = (0, — 0), 

which corresponds to a separation of ~ A/10. The tip has radius a = A/(107r), linear 
susceptibility = —0.4, second order susceptibility with ri^^^Eo/= 0.2 and all 

/ol 

other vanishing, and is scanned in the planes x = xq. In all numerical experiments, 
the incident electric held is a plane wave polarized in the x-direction with wave vector 
/c(r2)z. In Figure images of the extinguished power are shown in three diherent scan 
planes corresponding to xq = 2a, 2.5a and 3a. We see that the scatterers are well resolved 
with subwavelength separation in the closest scan plane and that the resolution degrades 
rapidly with distance from the plane x = 0. Qualitatively similar results are found for the 
intensity of scattered second harmonic light, as illustrated in Figure 
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Next we consider the case of THG. The setup is the same as in the case of SHG, except 
that the third order susceptibility of the tip is with = 0.2 and all other 

/qN 

Vijki vanishing. Once again, we hnd that the scatterers are better resolved in the closest 
scan plane. We also note that the relative extinguished power is smaller than in the case of 
SHG and that the intensity of the scattered third harmonic is correspondingly greater. See 
Figures!^ and 1^ In Figurethe extinguished power along a line in the closest scan plan is 
compared for SHG and THG. It is found that the the extinguished power for the case of a 
SHG tip is an order of magnitude larger than for a tip exhibiting THG. The corresponding 


result for the scattering amplitude is shown in Figure 10 


VII. DISCUSSION 

In this paper we have presented a generalization of the optical theorem for the scatter¬ 
ing of nonlinear electromagnetic waves. We consider in some detail the most important 
examples of quadratic and cubic nonlinearities. The theory is illustrated for the case of 
small inhomogeneities in the settings of second harmonic generation and the Kerr effect. 
In particular, we describe the manner in which scattering resonances are modihed by the 
presence of nonlinearity. As a second application, we consider the problem of computing the 
signal in a nonlinear near-held microscopy experiment. We note that the use of a tip with 
a nonlinear optical response affords the possibility of background suppression and spectral 
isolation of the detected signal. In future work, we plan to study the inverse problem for 
nonlinear near-held microscopy, whose goal is to reconstruct the linear optical properties of 
a sample illuminated by a nonlinear tip. This will necessitate the development of a scatter¬ 
ing theory that incorporates contributions from the tip and the sample and their respective 
interactions, as was done for the corresponding linear problem [26] . 


Appendix A: Conservation of Energy 

Here we show that for nondispersive media, the analog of the Manley-Rowe relations 
hold [8]. That is, if the susceptibilities Xi]\ xfji xifki purely real, then V ■ S = 0. 
We treat the cases of quadratic and cubic nonlinearity separately. Note that for an incident 


held consisting of a sum of a hnite number of frequencies, the integral in (28) becomes a 
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(a) Xq = A/(57r) 


(b) Xq = A/(47r) 


(c) xq = 3A/(IOtt) 


FIG. 5. (Color online) Apertureless NSOM images of the extinguished power for SHG. The 
susceptibility of the tip is = 0.2. Images are shown in the planes x = xq as indicated. 

The field of view of each image is 3A/(Svr) x 3A/(Svr). 



(a) Xq = A/(57r) (b) Xq = A/(47r) (c) Xq = 3A/(107r) 


FIG. 6. (Color online) Apertureless NSOM images of the extinguished power Pg for THG. The 
susceptibility of the tip is = 0.2. Images are shown in the planes x = xq as indicated. 

The field of view of each image is 3A/(Svr) x 3A/(Svr). 


sum. 








• • 


(a) Xq = A/(57r) (b) Xq = A/(47r) (c) Xq = 3A/(107r) 

FIG. 7. (Color online) Apertureless NSOM images of the far-held intensity at frequency 217 in 
SHG. The susceptibility of the tip is / 7 )^^^ = 0.2. Images are shown in the planes x = xq as 

indicated. The held of view of each image is 3A/(57r) x 3A/(57r). 



(a) Xq = A/(57r) (b) Xq = A/(47r) (c) Xq = 3A/(107r) 


FIG. 8. (Color online) Apertureless NSOM images of the far-held intensity at frequency 317 in 
THG. The susceptibility of the tip is = 0.2. Images are shown in the planes x = xq 

as indicated. The held of view of each image is 3A/(57r) x 3A/(57r). 


1. Quadratic Nonlinearity 


We begin by inserting the quadratic polarization (15) into the statement of energy con¬ 


servation (24). We then have 


V ■ S(r,a;) = 

+ f xifk(E(^i,(^2)E*(r,u)Ej(r,Ui)Ek(r,U2) 


(Al) 


^cji+a;2=G; 
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PJPq 



FIG. 9. (Color online) Extinguished power Pe/Po along the line defined by x = A/(57r) and z = 0, 
which corresponds to the closest scanning plane. Graphs are shown for SHG and THG. Here 
Po = a^cE'^- 


|Ap/|A„p 



FIG. 10. (Color online) Far-held intensity along the line dehned by x = A/(57r) and z = 0, which 
corresponds to the closest scanning plane. Graphs are shown for SHG and THG. Here Aq = apQ. 


19 






Making use of (28), the time-averaged divergence of energy current is 


^ • S(r) = {r]Uj)E*{r,u)Ej{r,u)) 

LO 

+ ^'^'^ix^ijli^]^^i,^2)E*{r,u)Ej{r,ui)Ek{r,U2)) . (A2) 

CJi+CJ2=^ 

The first sum is zero since Xif (i") ^)Ei (r, ^)Ej (r, uj) is real for each cn, due to the permuta¬ 
tion symmetry of The second sum vanishes as a consequence of the overall permutation 
( 2 ) 

symmetry of Xi/k constraint ooi + UJ 2 = oo- 

2. Cubic Nonlinearity 

For cubic nonlinearity we have 
^ • S(r)= ^^Im(xSj^(r;a;)E*(r,a;)Ej(r,a;)) 

UJ 

+ ^Iin(x!Ji(i';^i>^2,W3)E.(r,a;)Ej(r,a;i)Efc(r,a;2)Ez(r,a;3)) (A3) 

^1+^2+*^3=^^ 

The hrst sum is zero by the same argument as above. The second sum can be shown to 

( 3 ) 

be zero using the overall permutation symmetry of Xijki and the constraint wi -1-Ci ;2 + ti ;3 = 00 . 

Appendix B: Calculation of local fields 

Here we calculate the local fields of small scatterers for both second harmonic generation 
and the Kerr effect. 


1. Second harmonic generation 


We now calculate the local fields Ei{0,fl) and Ei(0,2fl). To proceed, we set r = 0 in 


(49) and (50) and thus obtain 


E,(0, Q) = fi) + k\Q)r]^]} / dVG,,(0, r'; Q)Ek{r\ Q) 

J |r^|<a 

+ 2A;2(f])pgj [ d^r'G,,{0, r'; fi)E,(r', 12)^;(r', 2fi) , (Bl) 

J |r^|<a 

E,(0, 2f2) = k^{2n)ri\l^ [ dVG,, (0, r'; 2fi)Efc(r', 2ff) 

J |r^|<a 

+ [ £r'G,jiO, r'; 2f2)Efc(r', f2)EKr', ff)) . 

J\r'\<a 


(B2) 


20 





Next, we use the fact that for a function gj 


* |r'|<a 


^ / / N / / N 47ra^ /I 2, , , 

<fr Gij{0,Y ]Uj)gj{Y ,uj) = —^ i- + i-k{uj) 


a‘^kP‘{oj) 


+ 0{k‘^{uj)a)] gi{0,uj) {1 + 0{k{uj)a)) , 

(B3) 


which is derived in Appendix B3 We then hnd that (Bl) and (B2) lead to a system of 


equations for the local helds which are of the form 

dvr 

li 


E,{0,n) = EincAOM) + ArkAn)AGR{n){gffEA0,n) + 2gl]lEA0,2n)El{0,n)) (B4) 


dvr. 


E,{0,2n) = ^kA2n)a^GR{2n){glfEA0,2n) + gA>EA0,n)Ek{0,n)) , 

where 


.( 2 ) 


(B5) 


(B6) 


The above is a set of nonlinear algebraic equations which we solve perturbatively. To 


proceed, we introduce a parameter e to scale the nonlinear terms in (B4) and (B5): 

E,{0,Q) = E,^,AO,Q) + ^kAQ)AGR{Q){vlfEA0,Q) + 2egi]lEA0,2Q)El{0,Q)){;B7) 


dvr. 


E,{0,2n) = ^kA2n)a^GR{2n){gffEA0,2n) + egl.lEA0,n)EA0,n)) ■ 

We then introduce asymptotic expansions for the helds of the form 

Bj(0,fi) = + (‘Ef\o,a) + ■ ■ ■ , 

£,(0,2fi) = £<“>(0, 2a) + e£i'>(0.2!J) + e"£p>(0.2f!) + • • • . 


,(2) 


(B8) 

(B9) 

(BIO) 


For simplicity, we consider the case of isotropic r/d) and 7 ]^) obeying permutation sym¬ 
metry. That is, gA = and = g^‘^\ with the other g^'^l vanishing. We also assume 

that the incident held Ej„c points in the a:-direction and the direction of observation s is 
taken to be in the ^(-direction. To simplify the notation, we set ffi = f2, f22 = 2ff, and write 


{Ei)j = Ej{0,kli), ki = Vti/c, {Einc)i = Einc,i(0,f2) and Gm = GAVti). Then and ( [B^ 
become 


(£,), = (£.„), + -^kla^Gm(ri'^\Ei), + 2er^g{El),(EM 

(E 2 ). = «(£,).+ £,(|(£.).(£,),) , (Bll) 
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Next, we expand the fields {Ei)j according to (B4) and (B5) and collect like powers of e. At 
0(1) we have that 


(£ 2 )!“’ = 0 . (B12) 


Thns 


mr = 


( 0 ) _ 


1 - 


At 0(e) we have 


dvr 


(B13) 


(Gj)!'’ = -kyGn^W'KE^f' + (£i)f) , (B14) 


which gives 


Thns 


(£ 2 )',“’ = ^kyGK2(r,^'HE,)i^^ + r,'SAEi)f\E,)f>) . 




At O(e^) we obtain 


I'p d 2 ) _ 2 3^ ,, ^ 

EfRiyij i^rz/i 


which gives 


Thus 


(BOf’ = y*:?a“G„,(r/‘>(Bi)f> + 2r,S(B2)S‘’(Br)f’) . 


,P 2(2) _ f *:p"G„i2r£>(B;)f>(B2)'i'' 

( 421)1 — - 


1 - f «:?o3Gri1,(» 


(B16) 


(B16) 


<S(B:)f*(B2)k) , (B17) 


(B18) 


(B19) 

We can now calculate the extinguished power Pe from ( |40[ ). We find that up to the order 

0(6^) 

8fli 


Pe = 


io=>Im ((R<‘'(Bi)i“> + R'‘'(£i)f> + 2R<t(£i-)S”’(£2)t’)(Bij;) . (B20) 
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2. Kerr effect 


Here we calculate the local field To proceed, we set r = 0 in (59) and thus 

obtain 


Ei{o, n) = EincA(^, n) + kAn)7]fJ / dA'CiAo, h; n)EAr^, n) 

J |r^|<a 

+ 34"(n)r,gL / <iVG.,(0.r'in)£;t(r'.!J)B,(r'.!2)B;(r',Sl)) . (B21) 

^|r'|<a 


We then find that (B21) leads to an equation for the local field which is of the form 


A'jr 

E,{0,n) =E,ncA^,^l) + —eaAkl)GR{p){A^^Ei{Q,^l) 
+?,il^lEAO,n)EAO,n)Ei{Q,n)) . 


(B22) 


The above is a nonlinear algebraic equation which we solve perturbatively. To proceed, we 


introduce a parameter e to scale the nonlinear terms in (B22): 


A'jr 

E,{Q,n) + —kA\Q)GR{Q){A^^Ei{0,n) 

+3e77g,E;(0,f2)Efc(0,n)E,(0,n)) . 


(B23) 


We then introduce asymptotic expansions for the field of the form 


Ei(0, n) = (0, n) + eEA’A, n) + e^E^’A, fi) + ■ ■ ■ 


(b/ 


,2 77'(2)/ 


(B24) 


For simplicity, we consider the case of isotropic and obeying the permutation 
symmetry. That is, r]A = and with all other rj^^h vanishing. We also 

assume that the incident field Ej„c points in the x-direction and the direction of observation 
s is taken to be in the 2 ;-direction. To simplify the notation, we write {E)j = Ej{0,Qi), 


k = Vl/c and Gr = Gr{Q). Then (B23) becomes 


dvr, 


(B), = (B.„)i + y+ 3e.,® (B*)j(B)t(B),) 


(B25) 


Next, we expand the fields {E)i according to (B24) and collect like powers of e. At 0(1) we 
have that 


(B)f> = + |^iV’G„r,“l(B)f' 


(B26) 
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Thus 


(E)f' = 


Ei, 


1 - ^k^a^GRT]W • 


At 0(e) we have 


(£)<■> = yfc^a=>G„(r,<‘)(B)|‘' + 3r,II,, 


which gives 


Thus 


and 




• 1 _ feWGR.,!') 




We can now calculate the extinguished power P^, from (40). We hnd that up 

0 (€) 


3. Evaluation of the Integral (B3) 


Here we show 


|r'|<a 


Atth^ /I 2 1 

d^r'Gij(0, r'; cj) = —-— ( - +i-k(cj) - ^ + 0(k(cj)E 


3 \a ' 3 ^ ' a^k‘^{u) 

For notational convenience, we put k = k{Q). Then 

(fr'Gij{0,r'-,uj) = [ d^r (Sij + ^didj \ -— 
^|r^|<a V ^ J T 


The second term vanishes when i j. Using the fact that 


Akr 


(V^ + -= —47r(5(r) , 


(B27) 

(B28) 

(B29) 

(B30) 

(B31) 
to the order 

:);) ■ (B32) 

(B33) 

(B34) 

(B35) 
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we have 


So 


Since 


r 1 pikr r 11/ 

/ dh-d,d,— = / d^r6,,-- -e— - 47rh(r) 

J\r\<a r Jirl<a ^ 3 \ r 


|r|<a 3r 3fc2y 


[ d^r^G,,{0,r'-oj)=6i,([ S 

'|r'|<a V>^|r|<a 


2 47r \ 


r-- 


r 


/ - = 

^|r|<a ^ 

we obtain the required result. 


|r|<a 3 r 3Fy 

47ra3 /I . 2 , 1 


7 + ^0^- -ir7w + 0(M , 


3 V a 3 a3fc2 


(B36) 


(B37) 


(B38) 


Appendix C: NSOM 

Here we derive the basic equations governing the NSOM experiments described in Sec¬ 
tion IVD 


1. Second Harmonic Generation 


The sample and the tip are taken to be small balls of radius a centered at fq, ri and r 2 . 
The corresponding susceptibilities are for k — ro| < a, = Vif for 

|r — ri| < a and |r — r 2 | < a, and xifki^'i^) = Vijl for |r| < o,- In this setting, the solutions 


to the wave equations of SHG (45) and (46) are 

Ei{r, O) = Einc,iir, O) -/ k^{Q)fi^]} 


'|r'-ro|<ci 


d^r'Gij{r,r';n)Ek{r',n) 


+ k‘^iSk)vfk J 

^ d^r'Gij{Y 

|r'—ri|<a 

+ ‘^k'^ik^Wjki 

[ d^r'G,j{ 

J |r'-ro|<a 


[ d?r'G,,{ 

7|r'-ro|<a 

+ k‘^{2n)r]f,} 

[ d?r'G,,{ 

J |r^—ri|<a 

+ k\2Q)7]f^i 

[ d^EG,,{ 

J\r'-ro\<a 




'\r'—r2\<a 


d?r'Gij{Y,Y']Vt)Ek{Y',Vt) 


(Cl) 


,(1) / 

J \r'-r 2 \<a 


d^r'GijiY, r'; C)Efc(r', 20) 


(C2) 
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Using the asymptotic form of the Green’s function given in (21), we hnd that the scattered 
helds are of the form 


E!{r,Q) = A{r,Q) 


^ik(Q)r 


r 

^ik{ 2 Q,)r 


(C3) 

(C4) 


uro 


Et{T,2n) = A{T,2n) 

where the scattering amplitudes are dehned by 

4(r, G) = - rA)k\n){fifjE,iTo, G) + 2G)E;(ro, 

+ - rA)k\n)vf^Ek{r,, , (C5) 

A(r,2G) = -f,f,)fc2(2G)(?}WE,(ro,2G) +77fjEfc(ro,G)Sz(ro,G))e*'’('')^^ 

(C6) 

Setting r = ro, r = ri and r = r2 in (Cl) and (C2), and carrying out the indicated 
integrations we obtain 


Ei{ro,Q) = ^*nc,i(ro,G) + -^k\Q)AGR{Q){4f E,{ro,Q) + 2r/gE,(ro, 2G)E*(ro, G)) , 

+ ^k^^iAa^Gijivi, ro; n)r]^]}Ek{ri, G) + ^k‘^{n)a^Gij{v2, fq; n)r]f^Ek{r2, G) (C7) 

E,(ro,2G) = yfc2(2G)a3GK(2G)(4')E,(ro,2G) +r7gE,(ro,G)^fc(ro,G)) 

+ yfc^(2G)a^Gy(ri,ro;2G)?7jJ,^Efc(ri,2G) + ^k‘^{2n)a^Gij{r2,ro-,2n)T]f^Ek{r2,2^) , 

(C8) 


Ei(ri,G) = Einc,i{ri,Q) + G) , 

+ YA:^(G)a^Gy (ri, fq; G)(i7jJ,^Efc(ro, G) + 2?7jJjEfc(ro, 2G)£;;(ro, G))) 

+ ^k\n)AGAri,r2;n)r]f-jEk{r2,n) (C9) 

U,(ri,2G) = ^k^{2n)AGR{2n)T]ffEJiT^,2^) 

+ ^fc^(2G)a^Gy (fi, fo; 2 kl){vfJEk{ro, 2G) + T]fJiEk{^o, kl)Ei{ro, G)) 

+ ^k\2n)AG,,iv,, V 2 ; 2n)rjfjEk{T2, 2G) , (CIO) 
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Ei{r2,Q) = Einc,i{r2,^) + -^k‘^{Q)a^GR{Q)r]l^'’Ej{r2,Q) , 

+ (rs, ro; n){fi\]}Ek{ro, Q) + 2r]^^J^Ek{ro, 2fi)£;;(ro, 12)) 

+ Y^^(^)«^<^*i(i’2,ri;12)77ji^Efc(ri,12) (Cll) 

E,(r2,212) = ^e{2Q)a^GH{2Q)vffE,{r2,2Q) 

+ ‘^k^{2Q)a^Gij{r2, fq ; 2Q){fi^^i}Ek{ro, 212 ) + r]ff}iEk{ro, Q)Ei{ro, 12 )) 

+ yA;2(212)a3G,,(r2,ri;212)4)Efc(ri,212) , (C12) 


Following the procedure indicated in Appendix B, the above equations can be solved per- 
turbatively for the local fields. 


2. Third Harmonic Generation 

As above, the sample and the tip are small balls of radius a centered at fq, fi and F 2 . 
The susceptibilities are x-j^(r;a;) = fj^j^ for |r — FqI < a, x^j^(r;a;) = for |r — ri| < a and 

/o\ ^o\ 

|f — r 2 | < a, and Xijkii^'i^) — VijLi kl — begin with the general cubic-nonlinear 

wave equations which are correct to order e: 

V X V X E(r,12) - fc2(12)E(r,12) =47rfc2(12)(x;j^(r, 12)E^(r, 12) 

+3xg,(r, 12,12, -12)E,(r, 12)E,(r, 12)E;(r, 12)) (.C13) 


V X V X E(r,312) -fc2(312)E(r,312) =47rA:2(312)(xSj^(r, 312)Ej(r, 312) 

+X%{r, 12,12,12)E,(r, 12)Efc(r, 12)EKr, 12)) (.C14) 


It follows immediately from (18) that the solution to (C13) and (C14) is given by 


E,(r,12) 


= E,„,,(r,12) + fc2(12) J dVxf,)(r',12)G,,(r,r';12)Efc(r',12) 

+ 3^2(12) J dV'xgL(r', 12,12, -12)Gi,(r, r'; 12)Efc(r', 12)EKr', 12)E;;(r', 12))(CT5) 
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Ei{r,3Q) = P{3Q) j 

+ 3k^3{Q.) [ ciVxgL(r', ^)G,,{r, r'; 3Q)Ek{r\ Q)Ei{r', Q)EUr\ f^)0Q16) 


For the specific set up described at the beginning of this section, (CIS) and (C16) become 


Ei{r, Q) = Einc,i{r, Q) + k‘^{Q)fifJ / (fr'Gij{r, r; Q)Ek{r, Q) 

J\r'\<a 


+ k‘^i^)vfk [ 

d^r'Gij{r,r'-,kl)Ek{r'A) 

J\r 

'-ri\<a 

+ [ 

d^r'Gij{r,r'-,kl)Ek{r'A) 

J\r 

'-r2\<a 

+ 3k‘^{Avfklm 

[ d^r'Gij{r,v']VL)Ek{r' 

J\r'\<a 

Ei(r,3fi) = k‘^{3A)vfj! j 

d^r'Gij{r, r'; 3r2)i7fc(r', 3Vt) 

r^|<a 

+ k‘^{3n)r]fj j 

d^r'Gij{v, r'; 3r2)i7fc(r', 312) 

r'—ri|<a 

+ k‘^{3n)r]f^ j 

d^r'Giji^r, r'; 3r2)i7fc(r', 312) 

r'—r 2 |<a 


(C17) 


+ k\3Q)vfL / d^r'G,j{r,r'-3Q)Ek{r',Q)Ei{E,Q)EUr^,Q)) . (C18) 

^|r'|<a 


Using the asymptotic form of the Green’s function given in (21), we find that the scattered 
held is of the form 

^ik(Q)r 


E!{r,n) = A{r,n)- 
E^{r,3A = 7li(r,3fi)‘ 


r 

^ik( 3 Q)r 


(C19) 

(C20) 


where the scattering amplitude is dehned by 


Ai(r, fl) = - fifj)k‘^{n){fj\l^EkiTo, n) + 3n^Efc(ro, ^)) 


47r 


jk 


Jk{fl)r-ri 


+ -^a (A -fifj)k (n)? 7 yEfc(ri,n)e 


(C21) 







Ai{r, 31]) = - rifj)k^{3Q){7]f^Ek{ro, 31]) + T]ff}i^Ek{ro, fl)Ei{ro, Q)Em{ro, 1])) 


jr-ri 


+ ^0=^(5,,■ -f,f,)fc2(31])4)E,(r2,31])e*'^(3n)f.r 


(C22) 


Setting r = ro, r = ri and r = r2 in (CIS) and (CIS), and carrying out the indicated 
integrations we find that 


Ei(ro,l]) = Einc,i{ro,Q) + -^a^k‘^{Q)filfGR{Q)Ej{ro,Q) 

+ Ya^A:^(l])?7ji^Gij(ro, ri; l])Efc(ri, 1]) 

+ ^a^k^{Q)ri^^jGij{r, r2; Q)Ek{r2, fl) 

+ ya3fc2(l])3r7;g,GH(l])i?,(ro,l])i?fc(ro,l])i?r(ro,l])) (C23) 

Ei{ri,n) = Einc,i{ri,n) + ^a^k^{Q)fi^l^Gij{ri,ro-,n)Ek{ro,n) 

+ ri; l])Ej(r', 1]) 

+ ^a^k‘^{Q)rifjGij{ri, r2; Q)Ek{r2, fl) 

+ Y^^^‘^^^">^'^EimGij{ri,ro;n)EkiTo,n)Ei{ro,n)E^{ro,n)) . (C24) 


Ei{r2,kl) = Einc,i{r2,kl) + —a^P{n)fiff}Gij{r2,ro;n)Ek{ro,n) 

+ ^a^fc^(l])?7]],^Gij(r2, ri; Q)Ek{ri, Q) 

+ ^a^k\n)rjffGR{n)E^iv2,n) 

+ ^a^k‘^{Q)3r]fJ^^Gij{r2,ro]fl)Ek{ro,Q)Ei{ro,Q)E^{ro,fl)) . (C25) 


E,{ro,3Q) = -^a^k\3Q)vlfGR{3Q)E,{ro,3Q) 

+ ri; 3Q)Ek{ri,3Q) 

+ ^a^k'^{3kl)r]fJGij{vo, r2; 3kl)Ek{r2, 31]) 

+ ^a^k‘^{3Q)r]fi}i^GR{3Q)Ek{ro,Q)Ei{ro,Q)Em{ro,kl)) ■ (C26) 
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Ei(ri,3fi) = —a^P{3n)fif^Gij{ri,ro-,3n)Ek{ro,3n) 

+ ^a^k\3n)r]ffGR{3n)E,{r,,3n) 

+ ^a^k^{3Q)r]^^f}Gij{ri,r2] 3Q)Ek{r2, 31]) 

+ ‘^a^k'^{3kl)rifJi^Gij{ri,ro-,3kl)Ek{ro,kl)Ei{ro,kl)Em{ro,kl)) . (C27) 

Ei(r2,31]) = Y«^^^(3^)^S^<^*i(i'2,ro;31])Efc(ro,31]) 

+ ‘^a^k'^{3kl)ri^^JGij{r2, ri; 31])Efc(ri, 31]) 

+ ^a^k\3n)vlfGR{3n)Ek{r2,3n) 

+ ^a^k^{3Q)r]fJi^Gij{r2, fq; 3Q)Ek{ro, Q)Ei{ro, Q)E^{ro, 1])) . (C28) 

Following the procedure indicated in Appendix B, the above equations can be solved per- 
turbatively for the local helds. 
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